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ABSTRACT 

Crack pr°p a ga tion tests were conducted to determine crack growth 

0 005 to o n n?n Ve h ® llc °P t ® r materials for surface cracks between 
0.005 to 0.020 inches m depth. Constant amplitude tests were 
conducted at stress ratios R = 0.1 and 0.5, and emphasis was 
placed on near threshold data (i.e. 10-8 to 10-6 inches/cycle). 
Spectrum tests were conducted using a helicopter spectrum. 

The test specimen was an unnotched tension specimen, and cracks 
were initiated from a small EDM notch. An opt ical /video system 
was used to monitor crack growth. The material for the test 
specimens was obtained from helicopter part forgings. Testincr 
was conducted at stresses below yield to reflect actual stresses 
m helicopter parts. 


I NTRODUC T I ON / B AC KG ROUND 


At the present time helicopters are generally "Safe Life" 
designed and maintained. This means that conservative structural 
component replacement times are defined based on SN fatigue data 
and parts are removed at these replacement times. Inspections 
may be performed at periodic intervals based on experience with 
service problems such as wear and corrosion. 

The Air Force damage tolerance approach, in which inspection 
intervals are quantified through crack growth analysis, has been 
applied to some helicopter structures to define inspection 
intervals in instances of service problems. Sikorsky Aircraft 
has also been investigating the general application of damage 
tolerance to the H-53 helicopter under Air Force contracts. 


* Test program performed by United Technologies Sikorsky 
Aircraft and Research Center under contract F09603-89-G- 
0096-0014. 
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As part of the Air Force damage tolerance investigations crack 
propagation data was generated for the actual materials used in 
the H-53 . Initial crack propagation data was obtained for large 
cracks using compact tension specimens, and for somewhat smaller 
corner cracks (greater than 0.015 inches) using edge notch 
specimens . 

This crack propagation data was then used to assess crack growth 
times and inspection intervals in H-53 rotor and airframe 
structure. These structures are subjected to significant 
vibratory loading at frequencies of the main and tail rotor rpm 
(approximately 3 and 12 hz, respectively), and sometimes at two 
to four times these frequencies . As a result crack growth 
analyses often involved initial crack depths at the limit of 
detectability (0.005 inches); and, even at these small crack 
sizes, crack growth times in some structures were relatively 
short . 

One concern in these crack growth assessments was the validity of 
the crack growth rate data at crack depths as small as 0.005 
inches. A. test program was therefore conducted under Air Force 
funding with the objective of obtaining crack growth rate data 
for surface cracks between 0.005 and 0.020 inches in depth in H- 
53 dynamic system and airframe materials. Emphasis was placed on 
constant amplitude near threshold data (i.e. 10-8 to 10-6 
inches/cycle). Some spectrum testing was also performed to 
examine load interaction effects. 


MATERIAL SELECTION AND DESCRIPTION 


Five materials were chosen for the test program based on a survey 
of materials in primary H-53 rotor and airframe structure. The 
materials are as follows: 

1. A17075-T73 forging 

2. 4340 (150 ksi) steel forging 

3. Ti-6A1-4V alpha-beta forged, annealed 

4. Ti-6A1-4V alpha-beta forged, solution heat treated and 
annealed 

5. Ti - 6A1 -4V beta STOA forged 

These materials encompass most of the H-53 structure, and are 
also common in many other Sikorsky helicopters such as the H-60 
series (Black Hawk, Sea Hawk, and Night Hawk) . The A17075-T73 
are used in a number of main and tail rotor parts, and 
^ost of the major frames and fittings in the airframe are 
fabricated from this material. The 4340 steel and three 
titaniums are used primarily in main and tail rotor parts. 

All the materials are die forgings and they were all obtained 
from actual helicopter parts bought to Sikorsky specifications as 
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verifi^t-hbf Conformance to the specifications was 

The usfof h artSi S tandard Sikorsky materials process evaluations 
essentL? £ in f ° r materlal sources was considered 

properties pr ° vldln9 pr °P er Serial form, microstructure, and 

Tensile ultimate and yield strengths for the materials are 

are V ”S» basis'- 2 - ^ A17075 - T7 3 and 4340 steel properties 

properties are properl e e s from Mil-Hdbk-5E. The titanium 

Sikorskv average of the properties measured in the 

biKorsKy material conformance tests. 


EXPERIMENTAL METHODS 


The test facility at the United Technologies Research Center 

includes a°ful3v S -,°t ® CF seriea "Interlaken" test system which 
includes a fully automated servo-hydraulic, closed loop test 

Series'^ 9 ™ od9l . 320 ° controller. A Universal Tes? Program 
control 3 ^!'lt lCh F driven software system, was used to 

loading th tSSt under both constant amplitude and spectrum 


Crack Detection And Monitoring 


An opt ical /video crack detection and monitoring system was 
selected and developed for this program. This method was 
selected over other methods such as replicate (refs. 1 2) and 
electric potential drop (ref. 6) based on efficiency in 
developing the test procedures and performing the tests. Similar 
optical methods for crack detection using a camera instead of a 
video system had been successfully used in references 3 , 4 , and 5 
The major components in this system consisted of a macrophoto 
lens, a CCD video camera, video recorder, strobed fiber optic 
lighting, and a video monitor. The use of a strobed lighting 
system provided the ability to capture the image of the crack at 
the maximum load without stopping the tests. It was therefore 
possible to video monitor and record the crack growth in real 
time. Crack lengths were measured on the video monitor which was 
ca ibrated against measurements taken in the microscope The 
crack length measurements were taken as the projected crack 
length along a straight line in the crack direction. 


Specimen Configuration 


In the selection of the test specimen the two types shown in 
Figure 2 (refs 1,2, 4, 5) were considered. Both of these specimens 
provided small test sections which facilitated the detection and 
tracking of naturally initiated cracks. As shown in Fig 2a the 
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notched specimen (refs. 1 , 2 ) consists of single edge notch nlat-p 

of radlufof d “ SiC T ° f d * 2 inches andla semfcirc^a/notch 
4 %? f ■ ' 125 lnC ^??‘ he unnotche d specimen in Figure 2b 
(rets. 45) is a smoothly tapered specimen, and therefore it has 
nosignrfrcant stress concentration. Its Overall dimensions a re 
6 ^ nches ; Although both specimens appeared reasonable the 

unnotched specimen was chosen for the following reasons: 

The smaller overall specimen size facilitated specimen 
manufacture from actual helicopter parts. 

It would not have to be stabilized at R = -1 testing as 
would the larger notched specimen. 

It had no stress concentration. This results in a 
simpler stress intensity solution. it was also 
believed that the notched specimen might require higher 
stresses at the crack origin to naturally initiate a 
crack . 


1 . 
2 . 
3 . 


The regions of the specimen in which cracks initiated were 

t ° mi " imize residual stresses and provide a surface 
suitable for visualization of small cracks. The specimen 
preparation involved an initial EDM process to form the taper 

withV? 4 n Sec ^ ion - Th ® specimens were then hand polished 
H r:^ 22 ° grit wet sand paper to remove the EDM remelt 

and heat affected areas. The next step was chemical milling of 
tne steel and. titanium specimens and electro polishing of the 

wlth^nn ^ e ?nnn ns ''l he £ inal P re P ara tion involved wet sanding 
with 6 °° to 1000 grit sand paper, and then polishing with 30 15 

and 3 micron diamond paste in a Varsol solution 


Test Method Development 


The test method development consisted of modifications to the 
specimen geometry, surface preparation, optical magnification, 
and specimen lighting. The specimen geometry was modified to 
that shown m Figure 3 to lower grip loads and thus prevent slip 
and failure at the grips. It was also determined that hiqhlv 
polished surfaces, and specialized fiber optic lighting were 
necessary for tracking the c rack. In early testing, naturally 
initiated cracks were attempted using a 30X magnification 
However, problems in detecting the crack led to a final 
magnification of 200x, and the use of a very small EDM starter 
notch to facilitate and localize crack initiation. Ver^ small 
EDM notches (0.006 to 0.009 inch surface lengths) were achieved 
by introducing an initial EDM notch, and then polishing the 
specimen to reduce the notch size. 
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EXPERIMENTAL PROGRAM 


Test Plan 


selected S.t? wa ? t0 P erform eight tests for each of the five 
selected H 53 materials. These tests consisted of two conqf^nt- 

0 mP ani U ^, teS ^ at each ° f three st ^ess levels R = a /a - -l 

constraints^he R° tests ' Due to budget and" technical 

constraints the R = -1 tests were not performed. The obiprHve 

0 005°to r 0 V 020 grow ^ h r ? te data for surface cracks from 

U.U05 to 0.020 inches m depth m the near threshold reaion Tt 

bes S t d poss?ble h Sata he but aCkS b t 1 -Rurally initiated to^&ain t£e 
vrZ t data, but, as noted, it was necessary to use small 

perform^he tests®? ‘ Aa ° ther ^Portant consideration was to 

levels tvDical o? str ? sses below yield to reflect stress 

levels cypical of helicopter structure. 


Stress Intensity Considerations 


One of the problems in the test 
propagation tests for 0 . 005 to 0 
intensity values. The problem i 
table shows the maximum stresses 
needed to initiate a crack in a 
to 1 million) ; the stress intens 
levels and crack depths of 0.005 
crack threshold stress intensiti 
intensities at the crack sizes o 
large crack thresholds. 


program was to perform 
.020 depth cracks at low stress 
s illustrated in Table 3. This 
(based on Sikorsky SN data) 
reasonable number of cycles (0.5 
ity ranges AK at these stress 
, 0.010, and 0.020; and the large 
es ; As shown the stress 
f interest are all well above the 


The method used to overcome this problem for stress ratios R = o 
Wa f "0. initiate the crack at R = -1 (compression fatigue) 
as suggested in reference 4. if the crack propagation tests are 

qtr^? er f° rmed A ^ R r 0 aad °' 5 without changing the maximum 
stress lower AK's than those in Table 3 may be realized. Using 
this method AK values near the large crack threshold were 

for a R n ^ d -l° r R = ° and °‘ 5 ‘ H °wever, no solution was obtained 


Another possibility for obtaining lower stress intensities is 
that the EDM starter notch might reduce the stress required f 
crack initiation from those shown in Table 3. However during 
development testing it was found that even with the EDM notch 
th S tbl CUlt t0 ^ nit:iate cracks at stresses much below those 


or 

the 

it 

in 
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Spectrum Considerations 


A number of helicopter stress spectra were examined for use in 

stanrtarrt h^' Cra ? k growth tests. These spectra included the two 
^f^ ard , hellC ° Pter rotor spectra HELIX28 and FELIX32 which were 
developed m a European collaborative effort for hinged and fixed 
rotor systems, respectively. Three spectra used in Sikorsky 
craclc propagation assessments of the H-53 were also examined. 
These spectra included those for the airframe pylon fold hinge 

s sis sis ie lbl * a - -<• 

Significant variability was found in these spectra. The HELIX32 
spectrum showed very little variation in maximum stress and its 
predominate stresses (stresses at which most cycles occur) are at 
j " ® FELIX28 spectrum has overloads in maximum stress 

around 1.5 times the predominate stresses, and its predominate 

overloadCr at R = 0 - 1 to 0 • 3 • The H-53 pylon fold spectrum has 
overloads in maximum stress of 1.2 to 2.4 times the predominate 

stresses, and its predominate stresses are at R values of 0.6 to 

Th * *~ 53 h ma l n rotor spindle has little variation in maximum 
e f, s ' an . lts stress ratios are mostly between 0.7 and 0.8. 

The H 53 mam rotor outboard blade spectrum has the most 
variability in maximum stress of all the spectra examined with 
some overloads as high as 3 to 4 times the predominate stress 

s P ectrum has a wide range of R values 
trom (J . 1 to 0.6 for high frequency stresses. 

In this test program a modified H-53 pylon fold spectrum was 
selected for use m the spectrum tests. The spectrum was 
modified to. maintain all R values below 0.6. This spectrum is 
representative of helicopter spectra with reasonably high 
overload factors . 


TEST DATA AND TEST RESULTS 


Table 4 provides a summary of the constant amplitude tests 
performed. At least two constant amplitude tests were performed 
tor each of the five materials at stress ratios R = 0.1 and 0 5 
The. cracks were initiated from the EDM notch in compression 
atigue (R =-1), and the maximum stress was either maintained or 
slightly increased for the crack propagation testing at R =0.1 or 
0.5. As explained, the initiation of cracks at R = -1 resulted 
in lower stress intensities than could be achieved by initiating 
the cracks at positive R ratios. The maximum stress is 
maintained the same for crack propagation to avoid overload 
effects . 

Table 5 provides a summary of the spectrum tests performed. One 
spectrum test was performed for each material using the H-53 
airframe pylon fold spectrum. The cracks were initiated using 
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constant amplitude compression fatigue (R = -1) . As shown in the 
table the maximum crack initiation stress was well below the 
maximum spectrum stress, and it was even below some of the high 
cycle maximum stresses in the spectrum. 

Crack length (a) versus cycle data (N) was recorded using the 
calibrated optical/video system for each constant amplitude and 
spectrum test. The constant amplitude data was converted to 
crack propagation rate data (da/dN vs Ak) using the ASTM secant 
method for determining da/dN, and a simple infinite plate 
solution for stress intensity K. The infinite plate solution 
assumes a semicircular crack, and examination of the fracture 
surface verified that this was a reasonable assumption. The 
mathematical expressions are: 


da/dN = (a i+1 - aj /N i+1 - N t ) 


and 


K = 0 . 7a (Tta) 1/2 

At the present time a limited amount of data analysis has been 
performed for the constant amplitude test. The results for R = 
0.1 are plotted in Figures 4 through 7. In all but Figure 4 
(A17075-T73) the small crack data is compared with large crack 
compact tension test data obtained about eight years ago in 
another Air Force funded test program. The A17075-T73 is 
compared with large crack propagation rate curves used in the 
crack propagation computer program FLAGRO. 


The AL7075-T73 R = 0.1 crack growth rate data (Fig 4) shows 
considerable scatter. The Al-2 specimen shows the most 
reasonable behavior with generally increasing da/dN with AK. The 
Al-16-1 and Al-3 specimens show erratic crack growth behavior. 

The Al-16-1 specimen generally showed a decreasing da/dN with 
increasing Ak. This may be due to the small initial crack size 
(Table 4) and interaction with the EDM notch. The Al-3-1 
specimen was tested at the highest Ak values and its behavior is 
also erratic. The Al-2 test specimen shows reasonable agreement 
with the FLAGRO curve. 

Some preliminary evaluation of the A17 07 5-T73 R = 0.5 data shows 
good agreement between the two specimens . The crack propagation 
rate da/dN is also generally increasing with increasing AK, but 
there is still a considerable amount of data scatter. 

The 4340 steel R =0.1 crack growth rate data (Fig 5) shows much 
less scatter than the aluminum, and remarkable agreement with the 
large crack compact tension threshold test results conducted 
eight years ago. The material sources for both the small and 
large crack specimens were H-53 main rotor hinge pins, but 
separate forgings were purchased from inventory and heat treated 
eight years apart . 
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Fiaure 6 I Sma11 CraCk R = 0 * 1 data is shown in 

6 a long with large crack compact tension threshold test 

^ he J ma11 crack data shows considerable scatter, and 
aster crack propagation rates than the large crack data. The 

' sour ^ es for th( r large and small crack specimens were a 
UH 60 mam rotor retention spindle forging, and a large H-53 
elastomeric mam rotor hub forging, respectively. 

Figure'7 A1 ;Lfwhh b f a R = sma11 c ^ack da ta is shown in 

lar ^ e crack compact tension threshold test 
data. As with the 4340 steel the Ti-6A1-4V alpha-beta crack 

r f e data /hows little scatter, and excellent agreement 
between large and small cracks. The small crack data contains 
both annealed and solution heat treated and annealed forginas 
and shows no difference in their crack growth rates. The large 

specimens were obtained from a large H-53 conventional main 

from^^smalle^^eteAt 9 ^ ^ Sma11 crack specimens were obtained 
irom smaller retention components (Table 1). 


CONCLUSIONS 


This test program has provided constant amplitude and snectmrr, 

iSches P in P deDth°? f ° r Sma11 surface cracks 0.005 to 0.020 

ncnes m depth for five common helicopter forged materials 

test spearmens were obtained from helicopter pfrt teg na and 

spIcT?icat?ons Wer ?he h te k t d a9ainSt Sikor ^ mSteriaHonlormance 
the yield stress to ^ Were conducted at stresses well below 
parts stress to be representative of stresses in helicopter 

indicates ^tha^the" 1 ^ 24 n ‘'l' 6 = ons ^ ant amplitude test data 

little crack arowfhrfto S H e t Tr-6A1-4V alpha-beta data show 

^ a crack growth rate data scatter when using the ASTM s pram 

method to calculate da/dN. These materials also showed aood 
?he e “75 W T73 and g ?' C ^? k 4, C ,T PaCt tenSlon thrlshoK^^ta. 

considerlb il Lf» d Z?‘ A1_4 Y beta ST0A materials show 

calculate dl/dN T When If ing the AS ™ secant method to 

nossiblv H,,o • J h aluminum data showed erratic behavior 
possibly due to interaction with the EDM notch Thp Ti fiai a\t 

s ss: as.'ss zss 

material^ Al r Force funded test program data for helicopter 

tensiondara ° th f . test programs include the compact 

corne? crack dataT" f ls paper ' and som e notched specimen 

corner crack data for initial crack radii between 0 015 and n mn 
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Further testing is recommended to obtain naturally initiated 
crack propagation data to avoid possible interaction with the EDM 
starter notch, and to perform tests at R = -1. It is also 
recommended that tests be performed at lower stress intensities 
to determine if 0.005 inch deep cracks will grow below the large 
crack threshold. 

Performing tests at lower stress intensities and at R = -1 for 

0.005 to 0.020 inch deep cracks is a difficult problem. A load 
shedding test was performed for one of the constant amplitude R = 
0.1 A1-7075-T73 specimens (specimen Al-16). In this test the 
loads were shed abruptly, and the data indicates that the crack 
propagated below the large crack threshold. The crack growth 
rates (da/dN) calculated by the secant method showed considerable 
data scatter. 
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Table 1. Materials For Small Cracks Program 


Material 


A17075 -T73 
Die Forging 


Steel 4340 (150 ksi) 

Die Forging 


Ti-6A1-4V Beta STOA 
Die Forging 

Ti-6A1-4V Alpha-Beta 
Die Forging , Annealed 

Ti - 6A1 - 4V Alpha -Beta 
Die Forging, Solution 
Heat Treated and 
Annealed 


Source 


SH-60 STA 295 
Side Fuselage 
Frame 

H-53 Main Rotor 
Horizontal Hinqe 
Pin 

H-53 Main Rotor 
Elastomeric Hub 

H-53 Main Rotor 
Hub 

S-76 Main Rotor 
Spindle/Cuff 


Part 

Numb er / Quan tity 

70209-222054-002 

1 


65103-11020-006 

2 


65103-11501-103 

1 

65102-11541-002 

2 

65102 - 11541-002 
1 


Table 2. Material Properties 


Material 

Ultimate 

Yield 


(ksi) 

(ksi) 

A17075 -T73 
Die Forging 

66 

54 

Steel 4340 (150 ksi) 

Die Forging 

150 

132 

Ti - 6A1 - 4V Beta STOA 
Die Forging 

133 

145 

Ti-6A1-4V Alpha-Beta 
Die Forging, Annealed 

139 

129 

Ti-6A1-4V Alpha -Beta 
Die Forging, Solution 
Heat Treated and 
Annealed 

155 

147 


i 



Table 3. Fatigue Stresses and Threshold Stress Intensity Factor 




Fatigue 


AK (ksi (Vin) 


Large Crack 


R 

Vibratory Stress 




AK t h 

Material 

Ratio 

(ksi) for 

a = .005 in 

a =.01 in 

a = .02 in 




Initiation in 




(ksi (Vin) 



0.5 to 1.0 x 10 6 





A17075-T73 

-1 

29-34 

5. 1-6.0 

7.2-8. 4 

10.1-11.9 


Forged 

0.0 

20-22 

3 5-3.9 

5. 0-5. 5 

7.0-7. 7 

2. 5-3. 5 


0.5 

13-14 

2. 3-2. 5 

3. 2-3. 5 

4.6-4. 9 

1.5-2. 2 

Steel 4340 

-1 

70-74 

12.3-13.0 

17.3-18.3 

24.6-26.0 


(150 ksi) 

0.0 

58-60 

10.2-10.5 

14.4-14.9 

20.4-21.0 

4-7 

Forged 

0.5 

35-36 

6 1-6.3 

8.7-89 

12.2-12.6 

3-4 

Ti-6A1-4V 

-1 

80-85 

14.0-14.9 

19.8-21.1 

28.0-29.8 


Beta STOA 

0 0 

55-59 

9 7-10.4 

13.6-14.6 

19.3-20.7 

4.6 

Forged 

0.5 

34-35 

6.0-6. 1 

8.4-8. 7 

1 1.9-12.3 

2. 5-4.0 

Ti-6A1-4V 

-1 

40-60 

7.0-10.5 

9.9-14.9 

14.0-21.0 


Alpha-Beta 

0.0 

30-50 

5. 3-8.7 

7 4-12.4 

10.5-17.5 

3-4 

Forged 

0.5 

20-30 

3. 5-5. 3 

5.0-7 4.2 

7.0-10.5 

2-3 


AK — 0.7(2Sv)V7ta 




Table 4 


Small Crack Constant Amplitude Test Summary 


Material 


A17075-T73 Forged 


Steel 4340 (150 ksi) Forged 


Ti-6A1~4V Beta STOA 
Forged 


Ti-6Al-4 Alpha-Beta Forged 
Annealed 


Ti-6A14V Alpha Beta 

Forged Solution 

Heat Treated & Annealed 


Specimen 

Number 

R 

Max 

Stress 

(ksi) 

Vibratory 

Stress 

(ksi) 

EDM 

Notch 

(in) 

Half Crack 
Length Range 
(c/2, inches) 

Delta K 
Range 

(ksi(in)**0.5) 

Al-3 

Al-16 

Al-2 

A1-R6 

A1-M7 

0.1 

0.1 

0.1 

0.5 

0.5 

45 

45 

30 

24 

24 

20.25 

20.25 

13.5 

6 

6 

0.0095 

0.005 

0.0075 

0.0085 

0.00875 

.0048-.0169 

.0025-.0056 

.0045.0229 

.0056-.-220 

•0048-.0217 

3.5- 6.S 

2. 5- 3. 8 
2.2-5.1 
1. 1-2.2 
1. 0-2.2 

SI -2 
SI -3 
Sl-1 

S3-1 

0.1 

0.1 

0.5 

0.5 

70 

70 

86 

66 

31.5 

31.5 

21.5 

16.5 

0.0077 

0.0068 

0.006 

0.0075 

.0042.0211 
.0047-.0217 
.0050-.0224 
.0061 -.02 18 

5.1-11.3 
5.3 11. 5 

3.0- 6.5 

3. 1- 6.0 


T 13 
Ti 4 
Ti 6 
T-20 
T 21 


1-1038-2A 

1-1013-1 

1-1038-5A 

1-1013-3A 


1-1037-1 
1-1037-4 
1-037-5 
1-037 9 


0.1 

0.1 

0.1 

0.5 

0.5 


0.1 

0.1 

0.5 

0.5 


0.1 

0.1 

0.5 

0.5 


80 

67.5 

70 

80 

75 


45 

50 

52.5 

50 


SO 

65 

60 

50 


36 

30.4 

31.5 
20 

18.75 


20.25 

22.5 
13.13 

12.5 


22.5 
29.25 
15 

12.5 


0.0093 

0.009 

0.0084 

0.0085 

0.0085 


0.0085 

0.00985 

0.00875 

0.008 


0.008 

0.0071 

0.008 

0.0085 


.0054.01 83 
.0053.02 18 
.0049.0229 
.0068.0238 
.0049-. 022 6 


.0075-.0225 
■0058-.0206 
.0075-.0205 
.0067.02 10 


.0068-.0228 
.0044.0230 
.0051-.0201 
.0065.02 10 


6.5- 12.1 

5.5- 11.1 
5.4 11.8 
4.1-7. 7 
3.3-7.0 


4.3-7. 5 
4.2-8.0 
2. 8-4. 7 
2. 5-4.5 


4.5- 8.4 
4.8-11.0 

2. 6- 5.3 
2.5-4.5 


Table 5. Small Crack Spectrum Test Summary 


Ug£UgitiatiQn—-- -> CrackPropagation 



Specimen 

Number 

EDM 

Notch 

(in) 

Max 

Stress 

(ksi) 

R 

Max 

Stress 

Half Crack 
Length Range 
(c/2, inches) 

Total 

Cycles 

A17075-T73 Forged 

AI-M6 

0.0078 

24 

-1 

53 

0044-.0204 

340179 

Steel 4340 (150 ksi) Forged 

Si -5 

0.0055 

70 

-1 

106 

0049-020 

5 1 5424 

Ti-6AJ-4V Alpha-Beta Forged 

T 1 - 1 2 

0.0085 

70 

-1 

106 

.0054-. 02 1 9 

53019 

Ti-61-4V /YJpha-Beta Forged Annealed 

1-1013-6 

0.0099 

55 

-1 

92 

0055- 0214 

220896 

Ti-6.\J-4Y Alpha-Beta Forged Solution Heat 
Treated and .Annealed 

1 -037-7 

0.00685 

57.5 

-1 

92 

.0049-. 0203 

533832 
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SUMMARY 


This paper presents recent developments in advanced analysis methods for the computation of stress 
intensity factors and the J-integral under conditions of small scale yielding in structural panels with multi- 
site damage. The method of solution is based on the p-version of the finite element method. Its implemen- 
tation was designed to permit extraction of linear stress intensity factors using a superconvergent extraction 
method (known as the contour integral method) and evaluation of the J-integral following an elastic-plastic 
analysis. Coarse meshes are adequate for obtaining accurate results supported by p-convergence data. The 
elastic-plastic analysis is based on the deformation theory of plasticity and the von Mises yield criterion. 

The model problem consists of an aluminum plate with six equally spaced holes and a crack emanating 
from each hole. The cracks are of different sizes. The panel is subjected to a remote tensile load. Experi- 
mental results are available for the panel. The plasticity analysis provided the same limit load as the exper- 
imentally determined load. The results of elastic-plastic analysis were compared with the results of linear 
elastic analysis in an effort to evaluate how plastic zone sizes influence the crack growth rates. The onset of 
net-section yielding was determined also. The results show that crack growth rate is accelerated by the 
presence of adjacent damage, and the critical crack size is shorter when the effects of plasticity are taken 
into consideration. This work also addresses the effects of alternative stress-strain laws: The elastic-ideafly- 
plastic material model is compared against the Ramberg-Osgood model. 


INTRODUCTION 

Reliable stress intensity factor computation of cracked structural components is of major importance in 
modern design of aircraft structures where requirements for residual strength and fatigue crack propagation 
must be met In this paper the problem of computing the stress intensity factors for a structural panel with 
multi-site damage, with guarantee of reliability, is discussed. This work focuses on the accuracy and reli- 
ability of the numerical solution of a proposed mathematical model and on the influence that different mod- 
eling assumptions may have on the results of the analysis. Two kinds of error are important in this case: 

1 . The differences between the exact solution of the mathematical problem formulated to represent a 
physical system or process and its numerical approximation are called errors of discretization. Is it 
possible to guarantee that the error of discretization Is small? 


1. Work supported by NASA Lyndon B. Johnson Space Center under Grant NAG 9-622. 


2. The differences between the exact solution of the mathematical problem formulated to represent a 
p ysical system and the actual response or behavior of the physical system are called errors of ideali- 
factors? 1 m ° dehng err0IS ' How does incol P° ra tmg plasticity in the analysis affect the stress intensity 

In some cases the two errors, the errors of discretization and modeling, may partially cancel one another. 
Therefore it is important to verify by means other than the experiment itself, that the numerical solution is 
close to the exact solution of the model. Only then is it possible to investigate whether the errors of ideali- 
zation are large or small by making comparisons with experimental observations. The efficient and reliable 
control of numerical errors, achieved by the use of a superconvergent method for the computation of the 

stress intensity factors, makes it feasible to investigate the sensitivity of crack extension to alternative mod- 
eling decisions. 

In this paper, the effect of plasticity on the values of the fracture mechanics parameters is investigated. 
Hrst, the finite element implementation is discussed; and second the method is illustrated with the integritv 
assessment of a row of fastener holes with multi-site damage, for which experimental results are available. 

FINITE ELEMENT IMPLEMENTATION 

Choice of an Extension Process 


In the finite element method the control of the errors of discretization can be achieved by mesh refine- 

men t (h- extensi °n) by increasing the polynomial degree of elements (p-extension), or a combination of 
both (hp-extension). 

The size of a finite element is the diameter of the smallest circle (or sphere) that contains the element. 
This diameter is denoted by h and the diameter of the largest finite element in the mesh is denoted bv h 
h-Extension involves letting h niax 0. Alternatively, we can hold the number of elements constant and*”** 
increase the polynomial degree of elements. The polynomial degree of elements is a vector p. p-Extension 
invo ves letting the smallest polynomial degree p min -+ °°. In hp-extensions mesh refinement is combined 
with an increase of the polynomial degree of elements. In these processes the number of equations that has 
to be solved, the number of degrees of freedom, is progressively increased, hence the name ‘extension*. 

Note that h- and p-extensions can be viewed as special cases of hp-extension, which is the general discreti- 
zation strategy of the finite element method. 

An important question is; “Which is the most efficient method of extension with respect to reducing dis- 
cretization errors? This question can be answered on the basis of a simple classification of the exact solu- 
tion of the problem one wishes to solve. The exact solutions (u EX ) have been classified into three main 
categories: 

Category A : u EX is analytic. A function is analytic at a point if it can be expanded into a Taylor series 
about that point on the entire solution domain, including the boundaries of the solution domain. Alterna- 
tively the domain can be divided into subdomains and u EX is analytic on each subdomain, including the 
boundary of each subdomain. The finite element mesh is so constructed that the boundaries of the subdo- 
mams are coincident with element boundaries. 

Category B: u EX is analytic on the entire domain, including the boundaries, with the exception of a finite 
number of points (in three dimensions u EX needs not be analytic along a finite number of lines). The 
mes is so constructed that the points where u EX is not analytic are nodal points (in three dimensions the 
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